The Power of Cosmic Flexion in Testing Modified Matter and Gravity 
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Flexion is the weak lensing effect responsible for the weakly skewed and arc-like appearance of 
lensed galaxies. The flexion signal-to-noise ratio can be an order of magnitude larger than that 
of shear. For the first time, we show how this makes flexion an invaluable tool for discriminating 
between alternative cosmological models. We analyse a scalar field model of unified dark matter and 
dark energy, a brane-world cosmology and two f(R) modified-action theories. We show that these 
models can be distinguished from ACDM at several standard deviations by measuring the power 
spectrum of cosmic flexion. 

PACS numbers: 98.80.-k, 98.80.Es, 95.36.+d, 95.36.+X 



Introduction. — In the last decades, cosmologists pro- 
posed several models alternative to the concordance A 
cold dark matter (ACDM) paradigm. These models at- 
tempt to find an agreement at least as good as that of 
ACDM with current cosmological datasets, such as the 
temperature anisotropy pattern of the cosmic microwave 
background radiation [1| , the dynamics of the large-scale 
structure of the Universe Q and the present-day cosmic 
accelerated expansion Q. However, in these theories, 
crucial topics such as the missing mass in galaxies and 
galaxy clusters and the current Universe's accelerated 
expansion are not explained by the usual dark matter 
and the cosmological constant A. On the contrary, these 
models mainly rely on either a modification of the law of 
gravity or the introduction of additional scalar or vector 
fields in the Universe's content. 

The family of alternative models with additional fields, 
also named modified matter models, includes dynam- 
ical dark energy or quintessence but also models 
which attempt to identify both the dark matter and 
dark energy effects with the properties of a single "dark 
fluid" ]|| ■ Conversely, the class of modified grav- 
ity includes a variety of approaches, which can how- 
ever be well represented by brane-world cosmologies 
Q and modified-action theories [||. Brane worlds de- 
scribe a four-dimensional "brane," which is our own Uni- 
verse, embedded into a higher-dimensional spacetime, the 
"bulk." In this scenario, Einstein's general relativity is 
still valid, but the higher-dimensional behaviour of grav- 
ity induces non-negligible signatures on the Universe's 
evolution and the growth of cosmic structures on the 
brane. Finally, modified-action theories directly modify 
the law of gravity by generalising the Hilbert-Einstein La- 
grangian. Among all the possible theories, f(R) gravity, 
where Ricci's scalar R is replaced by a generic function 
f(R), is probably the most investigated approach. 

In this Letter, we choose three models to explore the 



space of modified matter and gravity theories. Specif- 
ically, we consider a model of unified dark matter and 
dark energy Q , a phenomenological extension of the well- 
known DGP brane-world cosmology 0, and twoJ(R) 
models able to pass the Solar system gravity tests [lll . ll2| . 
All of them reproduce the ACDM expansion history, thus 
representing viable alternatives for the description of the 
background evolution of the Universe. To be able to dis- 
criminate between them and ACDM it is therefore crucial 
to investigate the regime of cosmological perturbations. 
This analysis has been carried out using several observ- 
ables, for instance the power spectrum of density fluc- 
tuations and cosmic shear However, it is not rare 
that the predicted signal is very similar to what expected 
in ACDM. Here, we show that the degeneracy between 
models can be lifted by cosmic flexion, namely the flex- 
ion correlation function whose signal originates from the 
large-scale structure of the Universe. We will present 
parameter forecasts and additional gravitational lensing 
statistics elsewhere 
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Cosmic Flexion. — The deflecting gravitational field of 
the extended large-scale structure of the Universe - which 
is simply the Newtonian potential $, in general relativ- 
ity - is responsible for deflection of light rays emitted 
by distant sources. This phenomenon is known as weak 
gravitational lensing. Therefore, photon paths from a 
galaxy located at 9 on the sky are deflected by an angle 



a 



dip, 



(1) 



where d — d\ + id% is the gradient with respect to di- 
rections perpendicular to the line of sight and ip is the 
projected deflecting potential. Unfortunately, the deflec- 
tion angle is not observable directly. This is because one 
does not know the true two-dimensional distribution of 
the sources on the sky. On the other hand, its gradient, 
the distortion matrix d a dbi>, is measurable. In particular, 
the entries of the distortion matrix can be related to the 
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effects of convergence k and (compfex) shear 7 occurring 
to the source image, i.e. 



k = ^dd*ip, 7 = hs&tj). 



(2) 



If convergence and shear are effectively constant within 
a source galaxy image, the galaxy transformation is 
0' a = A ab 9 b , where A ab = 6 ab -d a d b ip and a,b = 1, 2 label 
the coordinates on the sky. Flexion arises from the fact 
that the shear and convergence are actually not constant 
within the image, it therefore represents local variabil- 
ity in the shear field that expresses itself as second-order 
distortions in the coordinate transformation between un- 
lensed and lensed images. Thus, by expanding the ob- 
served galaxy position 9' at the second-order in the de- 
flection angle, it follows that 9' a — A ab 9 b + D abc 9 b 9 c /2 
15], with D abc = d c A ab . As the distortion matrix can 
be decomposed into the convergence and the shear, it is 
usual to define a spin-1 and a spin-3 flexion, which read 



J = -ddd*ip, S = -dddtp, 



(3) 



resp ectively. Since measurements of S are noisier than IF 
[161 ] , we will restrict our analysis to 5" only. 

To construct the flexion correlation function from 
large-scale structure, we start from the definition of the 
projected deflecting potential, 



<l>(0) = Jd x w(xMx,o), 



(4) 



where dx = dz/H(z) is the radial comoving distance, 
H(z) is the expansion history of the Universe and W(x) 
is the weak lensing selection function 17|. W(\) depends 
on the redshift distribution of the sources n[x(z)], nor- 
malised such that fdxn(x) = 1- 

In the flat-sky approximation, we expand the flexion 
in its Fourier modes Hence, from the definition of 

angular power spectrum 



{?{l)T (£')) = (2ir) 2 S D (£ - t')C s {t), 



(5) 



which is the Fourier transform of the two-dimensional 
correlation function, and from Eq. we finally get 18] 



(6) 



Modified Matter/ Gravity Models. — We now briefly re- 
view the three models we use. We refer to them as: UDM 
for the model of unified dark matter and dark energy; 
eDGP for the phenomenologically extended DGP brane 
world; and St and HS for the two f(R) theories of [ll[ 
and [12]. respectively. 

In the class of UDM models we use H, a single scalar 
field with a Born-Infeld kinetic term [19j mimics both 
dark matter and dark energy. The energy density of the 



scalar field reads pudm = Pbm + PA, where Pdm oc a~ 3 
and pa = const., that yields the ACDM Hubble pa- 
rameter exactly. However, there also is a pressure term 
Pudm = — PA which leads to a non-negligible speed of 
sound for the perturbations of the scalar field itself. This 
is a common feature in modified matter models and it 
typically causes an integrated Sachs- Wolfe effect incom- 
patible with current observations [20( . To solve this prob- 
lem, we use the technique outlined in Q, where the au- 
thors construct a UDM model able to reproduce both 
the correct temperature power spectrum of the cosmic 
microwave background and the clustering properties of 
the large-scale structure we see today. 

The sound speed is parameterised by its late-time value 
Cqo (in units of c = 1) and the growth of cosmic struc- 
tures strongly depends on it. Indeed, the presence of the 
sound speed produces an effective Jeans length Aj for 
the Newtonian potential. Thus, its evolution is no more 
scale independent. Specifically, the Fourier modes $fc 
are suppressed on scales k > 1/Xj and oscillate around 
zero. The larger is the value of Coo, the earlier the Newto- 
nian potential starts decreasing (for a fixed scale) or at a 
greater scale (for a fixed epoch) [21| . Since the Newtonian 
potential is responsible for light deflection, weak lensing 
is a powerful tool to constrain UDM models 2l| and in 
particular three-dimensional cosmic shear [22j |. However, 
UDM models with Coo < 10~ 3 still produce a signal vir- 
tually indistinguishable from that of ACDM. 

In the eDGP model [lOj, the cross-over length r c , which 
defines the scale at which higher-dimensional gravita- 
tional effects become important, is tuned by a free pa- 
rameter a G [0,1]. It is strictly related to the graviton 
propagator. Particularly, a = and a = 1/2 reduce to 
ACDM and standard DGP, respectively. Recently, it has 
been shown that the eDGP model excellently fits geo- 
metrical datasets such as the Hubble diagram of type la 
supcrnovae and gamma ray bursts, the scale of baryon 
acoustic oscillations and the CMB distance indicators 
,23|. Therefore, it is crucial to test this model in the 
regime of cosmological perturbations. 

As in other modified gravity theories, the two metric 
perturbations, the Newtonian potential $ and the metric 
potential evolve differently even with no anisotropic 
stress. Contrarily, for general relativity $ = — ^ holds in 
the matter dominated era. Thus, when we study gravi- 
tational lensing we have to deal with the deflecting po- 
tential T = (fy — $)/2. Moreover, its Poisson equation, 
which relates it to the distribution of the matter overden- 
sities, is modified by the presence of an effective time- and 
scale-dependent gravitational constant. 

It is worth giving a final remark on the evolution of 
matter fluctuations. Unlike the linear growth of pertur- 
bations, that can be described analytically, the non-linear 
regime has to be explored numerically. Two approaches 
have been followed and we refer to them as "KW" and 
"PPF." The former generalises the halofit procedure 
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24j to the eDGP scenario according to recently per- 
formed iV-body simulations 25[, whilst the latter in- 
terpolates the eDGP non-linear matter power spectrum 
with that of ACDM in order to reproduce general rel- 
ativity at small scales and be thus able to pass Solar 
system gravity tests j26[. The functional forms of this 
last approach have been obtained by perturbation theory 
27] and confirmed by TV-body simulations [28]. Unfortu- 
nately, were PPF the correct non-linear prescription, we 
would not be able to discriminate between the eDGP and 
ACDM signals even with the present and next generation 
weak lensing surveys (23[. 

Finally, we analyse the St and HS f(R) theories, which 
also are degenerate with ACDM at background level . 
Their functional form allow them to achieve the late-time 
accelerated expansion of the Universe with no formal cos- 
mological constant. On the contrary, they present three 
free parameters, c±, C2 and n. It has been shown that 
the growth of linear perturbations strongly depends on 
the function f(R), which acts by generating a time- and 
scale-dependent gravitational constant, as well as an ef- 
fective anisotropic stress (30| . Regarding the non- linear 
evolution of perturbations, the PPF technique is still 
valid, as confirmed by A^-boy simulations [3l|. Cosmic 
shear studies on these models gave interesting results 
[32I [33I ] . but the St signal is nonetheless almost com- 
pletely degenerate with ACDM [34 1 . 

Results and Discussion. — Here, we present the cos- 
mic flexion power spectrum ([6]) expected in the al- 
ternative models outlined above and we compare it 
with the ACDM prediction. For this, we use a fidu- 
cial flat Universe where the Hubble constant is Hq = 
100 h km s _1 Mpc -1 and h — 0.7. The matter density in 
units of the critical density is fl m = 51dm + ^6 = 0.28, 
with f^DM and ilb — 2.22 • 10~ 2 h~ 2 the dark matter 
and baryon fractions, respectively. The tilt of the pri- 
mordial matter power spectrum is n s = 0.96 and the 
density fluctuation mis on the scale of 8 ft. -1 Mpc is 
a$ = 0.8. For the UDM model, we probe Coo = 5 • 10~ 4 
and Cqo = 10~ 3 . The eDGP model parameters are 
a = 0.116 and r c H = 155.041 [23j]. The St(HS) parame- 
ters read log 10 ci = 2.38(4.98), log 10 c 2 = -2.6(3.79) and 
n = 1.79(1.67) 0. 

It is important to note that there currently is no linear- 
to-non-linear mapping in UDM models. Nevertheless, 
differences between ACDM and UDM models arise at 
scales smaller than the sound horizon. With a cross-over 
wavenumber k ~ 1/Aj, if the sound speed is small enough 
to guarantee that Aj is well within the non-linear regime, 
we can assume that the non-linear evolution of the UDM 



power spectrum will be similar to that of ACDM [22 1 . 

We use the specifics of the upcoming ES A Euclid satel- 
lite faiHI Euclid is one of the ESA Cosmic Vision 2010- 



2015 approved projects and is currently in the timeline of 
M-class missions. Its survey area will be 20, 000 square 
degree, with a sky coverage / s k y — 0.48 and a source 
distribution over redshifts 36] 



n(z) oc z 2 e 



(7) 



where zq = z m /lA and z m = 0.9 is the median redshift 
of the survey. The number density of the sources, with 
redshift and shape estimates, is n = 35arcmin -2 . To 
compute errorbars, we use 



(2e + i)f sky 



[C*W+Nf\, (8) 



generalising thus the approach of [17J. This is because 
unlike shear - flexion has a dimension of length -1 (or 
ang kT 1 ). This means that the effect by flexion depends 
on the source size. Recently, it has been shown that the 
noise power spectrum Nf for flexion is inversely propor- 
tional to the squared angular scale [37} ; we therefore set 



Nl " l 2 n ' 



(9) 



with (5i n t 2 } ' 5 — 0.03arcsec 1 the galaxy-intrinsic flex- 
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1 http : //sci . esa. int/science-e/www/area/index. cf m?f areaid=102 



FIG. 1. Power spectrum of cosmic flexion C J (£) for the al- 
ternative cosmological models presented in the text. 

Fig. [1] shows the cosmic flexion power spectra 0^(1) 
of ACDM (red, solid), eDGP (green) with both KW 
(dashed) and PPF (dot-dashed) linear-to-non-linear 
mappings and the f(R) models (red) of St (dashed) and 
HS (dot-dashed). As expected, the UDM signal is sup- 
pressed at small angular scales because of the presence of 
the scalar field sound speed. The eDGP model is still very 
close to the ACDM prediction, particularly the PPF non- 
linear power spectrum, for it being specifically designed 
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to reproduce general relativity on small scales. On the 
other hand, the St and HS models clearly show the scale 
dependence of the Newtonian gravitational constant G. 
Indeed, in the so-called "scalar-tensor" regime it reaches 
the value - 4G/3 0. 

Nevertheless, cosmic flexion shows an outstanding im- 
provement in the separation between the signals when 
compared to the cosmic shear power spectrum. Indeed, 
the dark-grey shaded area represents the lcr-error region, 
whilst light-grey refers to errors six times larger. Flex- 
ion measurements are made on the shapes of the source 
galaxies, exactly as in the cosmic shear analysis. There- 
fore, the source number density n is the same for the 
two observables and, with a space-based, wide-field sur- 
vey such as Euclid, we can collect a fairly large statistics. 
However, the intrinsic flexion rms (^int 2 ) 0,5 is an order 
of magnitude smaller than the cosmic shear rms and the 
power spectrum is thus significantly less noisy. 

We conclude that cosmic flexion is an excellent tool for 
testing alternative cosmological models and discriminate 
between them. With realistic values for the mean galaxy 
number density n and the flexion noise Nf , which in- 
cludes its angular scale-dependence [37j . expected from 
the upcoming Euclid mission, we find an admirable sep- 
aration between cosmic flexion power spectra C 3 ^) of 
viable models which are almost degenerate with ACDM 
when investigated with other observables, such as cosmic 
shear. We will provide a more detailed analysis of these 
outstanding results in 14| . 
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